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ABSTRACT 

The relative strength between forward and reverse shock emission in early 
gamma-ray burst afterglow reflects that of magnetic energy densities in the two 
shock regions. We numerically show that with the current standard treatment, 
the fireball magnetization is underestimated by up to two orders of magnitude. 
This discrepancy is especially large in the sub-relativistic reverse shock regime 
(i.e. the thin shell and intermediate regime) where most optical flashes were 
detected. We provide new analytic estimates of the reverse shock emission based 
on a better shock approximation, which well describe numerical results in the 
intermediate regime. We show that the reverse shock temperature at the onset 
of afterglow is constant, (f^ — 1) ~ 8 x 1CT 2 , when the dimensionless parameter 
£o is more than several. Our approach is applied to case studies of GRB 990123 
and 090102, and we find that magnetic fields in the fireballs are even stronger 
than previously believed. 

Subject headings: gamma rays: bursts — magnetic fields — hydrodynamics — 
ISM: jets and outflows 



1. INTRODUCTION 

A widely accepted model for producing gamma-ray bursts (GRBs) is based on the dis- 
sipation of a relativistic outflow (e.g. Piran 2004; Zhang & Meszaros 2004). The internal 
energy produced by shocks is believed to be radiated via synchrotron emission. Although 
the presence of strong magnetic fields is crucial in the model, the origin and its role in the 
dynamics are still unknown. Understanding the nature of the relativistic outflow, especially 
the energy content, acceleration and collimation, is a major focus of international theoretical 



Astrophysics Research Institute, Liverpool John Moores University, Twelve Quays House, Birkenhead, 
CH41 1LD, UK 



*rmh@astro. livjm.ac.uk 



-2 - 



and observational efforts. Relativistic outflow from a GRB central engine is convention- 
ally assumed to be a baryonic jet, producing synchrotron emission with tangled magnetic 
fields generated locally by instabilities in shocks (Medvedev & Loeb 1999; Nishikawa et al. 
2005; Spitkovsky 2008). Recently an alternative magnetic model is attracting attention from 
researchers (e.g. Drenkhahn & Spruit 2002; Zhang & Kobayashi 2005; Lyutikov 2006; Gi- 
annios 2008; Mimica et al. 2009, 2010; Zhang & Pe'er 2009; Zhang & Yan 2011; Narayan et 
al. 2011; Granot 2012). The rotation of a black hole and an accretion disk might cause a 
helical outgoing Magnetohydrodynamic (MHD) wave which accelerates material frozen into 
the field lines (Tchekhovskoy et al. 2008; McKinney & Blandford 2009; Komissarov et al. 
2009). In the magnetic model, a fireball is expected to be endowed with primordial magnetic 
fields from the the central engine. 

The first detection of ten percent polarization of an optical afterglow just 160 sec after 
the GRB explosion (Steele et al. 2009) opens the exciting possibility of directly measuring 
the magnetic properties of the GRB outflow. Recently polarization measurements of the 
prompt gamma-ray emission were also reported (Kalemci et al. 2007; McGlynn et al. 2007; 
Gotz et al. 2009; Yonetoku et al. 2011). Although these polarization measurements suggest 
that at least some GRB outflows contain ordered magnetic fields and they are still baryonic, 
the sample is small and further observations will be necessary to confirm the magnetic model 
and/or to understand the role of magnetic fields in the dynamics. In this paper, we revisit 
the magnetization estimate of the GRB outflow (hereafter "fireball") based on photometric 
observations of early optical afterglow. It is more sensitive to the magnetic energy density, 
rather than the length scale of magnetic fields in the fireball, and it is complementary to 
polarimetric methods (e.g. Lazzati 2006; Toma et al. 2009). 

A steep decay in early optical afterglow light curves is usually considered as a signature 
of the reverse shock emission (e.g. Akerlof et al.1999; Sari & Piran 1999; Meszaros & Rees 
1999; Soderberg & Ramirez- Ruiz 2002; Li et al. 2003; Fox et al. 2003; Nakar & Piran 2005). 
The early emission contains precious information on the original ejecta from the central 
engine. The magnetization of the fireball can be evaluated by using the relative strength 
of the forward and reverse shock emission (Fan et al. 2002; Zhang et al. 2003; Kumar & 
Panaitescu 2003; Gomboc et al. 2008). However, the current method is obtained by using a 
simplified shock dynamics model which is known to be inaccurate in the intermediate regime 
between the thin and thick shell extremes (Nakar & Piran 2004), and actually most observed 
events occurred in the intermediate regime. Here we numerically re-examine the interplay 
between the forward and reverse shock emission at the onset of afterglow. In Section [2] we 
set out a simple conventional approach to understanding the two shock emissions and refine 
the definition of the magnetization parameter in Section [31 In Section H] we consider a new 
approximation to discuss the reverse shock emission in the intermediate regime. In Section 
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|5l we then test these analytic approximations with numerical simulations. In Section [6] we 
present case studies of GRB 990123 and 090102 in terms of the magnetization parameter. 
Finally in Section [7] we summarize the results. 



2. Forward and Reverse Shock 

We consider a homogeneous fireball of energy E and a baryonic load of total mass 
M confined initially in a sphere of radius r$. We define the dimensionless entropy To = 
E/Mc 2 ^> 1. This fireball expands into a homogeneous interstellar medium (ISM) of particle 
density n\. This can be considered to be a free expansion in its initial stage. After a short 
acceleration phase, the motion becomes highly relativistic, a narrow shell is formed. After 
the fireball shell uses up all its internal energy, it coasts with a Lorentz factor of To and the 
radial width A ~ r . 

The deceleration process of the shell is described with two shocks: a forward shock 
propagating into the ISM and a reverse shock propagating into the shell. The forward 
shock is always ultra relativistic, while the evolution of the reverse shock is determined by 
a dimensionless parameter £ = (1/ Ao) 1 / 2 ^^ 3 where I = (?>E / Anrripnic 2 ) 1 ^ is the Sedov 
length and m p is the proton mass. If £o < 1 ( so called thick shell case), the reverse shock 
becomes relativistic in the frame of the unshocked shell material and it drastically decelerates 
the shell. If £o > 1, the reverse shock is inefficient at slowing down the shell. The deceleration 
radius r d and the Lorentz factor of the shocked material at r d are usually approximated 
as r d ~ / 3 / 4 Aq /4 and T d ~ (Z/A ) 3 / 8 for £0 < 1, and r d ~ l/T 2 /3 and T d ~ T for £ > 1 (Sari 
and Piran 1995; Kobayashi et al. 1999). After the deceleration, the profile of the shocked 
ISM medium begins to approach the Blandford Sz McKee (1977) solution. 

We first discuss the forward and reverse shock emission by using these conventional 
estimates, and the accuracy (i.e. correction factors) will be numerically examined later. The 
deceleration of a shell happens at an observer time 

CI Q 

where Ct ~ (T d /To)~ 8 / 3 ~ max(l,£^" 2 ) and all the correction factors in this paper, including 
C t , are defined as ones relative to the conventional thin shell estimates. At the deceleration 



1 The details of the acceleration mechanism is not important in this paper, the same discussion can be 
applied even if a non-thermal process is responsible for the production of a relativistically expanding shell. 
Because of the relativistic beaming effect, the radiation from a jet before the jet break can be described by 
a spherical model. 
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time, the forward and reverse shock regions have almost the same Lorentz factor and internal 
energy density. However, the reverse shock region has a much larger mass density and 
therefore it has a lower temperature. Introducing a magnetization parameter R B = tB,r/ e Bj, 
it is shown that the typical frequencies v m and peak fluxes F vmax of the synchrotron emissions 
from the two shocks are related as (Kobayashi & Zhang 2003; Zhang et al. 2003) 

^m,r(^d) _ /-y r -2pl/2 Fv,max,r{td) r< r D 1 / 2 (o\ 

7TT — °mJ- o B > — p — 0-^8 > l Z J 

where C m ~ (r^/To) -4 ~ max(l,£^" 3 ) and Cf ~ (I^/To) 2 ~ min(l, Co^ 2 ) are correction 
factors, the subscripts r and / indicate reverse and forward shock, respectively. We have 
assume that the electron equipartition parameter e e and the electron power-law index p are 
the same for the two shock regions, but with different magnetic equipartition parameter 
tB as parametrized by Rg. The reason we introduce the Rb parameter is that the fireball 
might be endowed with primordial magnetic fields from the the central engine. We can give 
a simple relation u c>r / v c j ~ Rb^^ between the cooling break frequencies of the two shock 
emissions (Zhang et al. 2003). As we will see in the next section, this simple estimate is 
good enough for the magnetization estimate. 



3. Magnetization Estimate 

For no (or moderate) primordial magnetic fields in a fireball, we expect v m ^ r <C u m j and 
v c ,r ~ Vc,f at the peak time of the reverse shock emission t p ~ td- The optical band v ovt should 
satisfy a relation v m>r < v apt < v c>r during the early steep decay phase of the reverse shock 
emission. Otherwise the decay is much slower or faster than the typical decay t~ 2 (Kobayashi 
2000). There are four possible relations between the break frequencies and the optical band 
at the peak time t p : (a) v m>r < u opt < v m j < v CjT < v c j, (b) u m , r < v opt < v CjT < u mJ < u c j, 
(c) u mjr < u opt < y c ,r < v c j < v m j, and (d) v m , r < u mJ < v opt < v c>r < v c j. In the case 
(a) and (b), the forward shock emission peaks at t p j when the typical frequency u m j goes 
through the optical band (the top panel in figured]). Using v m j oc t~ 3 / 2 , we get the peak 
time and peak flux ratio 

Rt = tp,f/tp = { l 'm,}(tp)/ l 'a P t) 2Z i (3) 
Rf = FpJ Fp,f (Fv,max,r(tp) / Fu,max,f)^opt / ^m^itp)) ^ ^ (4) 

where F p and F p j are optical peaks in the time domain, while F V)maXjr and F vmax j are peaks 
in the spectral domain for a given time. The hydrodynamics evolution of a reverse shocked 
shell is investigated in Kobayashi & Sari (2000), and the decay index a ~ (3p + l)/4 ~ 2 of 
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the reverse shock emission is found to be almost independent of £o when u m>r < v opt < u, 
Combining Equations (J2J), © and (j4]), we obtain (Gomboc et al. 2008) 



c,r- 



D3 F 4a-7 \ 2 /(2a+l) 2/5 



At this stage, we assume that T is a known quantity, and we will discuss how to estimate 
T from early afterglow observations in section 5.3. In the case (c), if the forward shock 
emission makes a transition from the fast cooling to the slow cooling regime before it peaks, 
it becomes equivalent to the case (b). The estimate (JSJ) is valid. On the other hand, if it is 
still in the fast cooling regime when u c j crosses the optical band, the forward shock emission 
rises and decays very slowly as t 1//6 and t -1 / 4 (Sari et al. 1998). Since this behavior is not 
consistent with most early afterglow, we do not discuss the details q Finally, in the case (d), 
the forward shock emission also peaks at the onset of afterglow, and it follows that Rt — 1. 
It is possible to show that Equation ([5]) is still valid. 

When an early afterglow light curve shows a flattening at t = tfi at after the steep decay 
phase (the bottom panel in figure [1]), the reverse shock emission dominates at early times. 
The forward shock peak is masked by the reverse shock emission, the peak (t p j, F p j) is not 
observationally determined. In such a case, the upper limit t p j = tfi at gives a rough estimate 
of Rb- Considering that the reverse and forward shock emission components are comparable 
at the flattening, we obtain R F ~ R". Substituting this relation into Equation §5§, we get 

where R t = tfi at /t p . If the forward shock emission peaks earlier t p j < tfi at , the real value of 
Rb might be slightly different. To evaluate how Rb depends on t p j, we refer to the scalings 
R t oc t p j and Rf oc where «/ is the decay index of the forward shock emission. Using 
these scaling, one finds that the dependence is weak: Rb oc ^ +af a )/( 1+2a )_ jf f orwar d 
shock decays as the theory suggests a/ = 3(p — l)/4, a relation — a = holds, and 

Rb does not depend on t p j (Gomboc et al. 2008). 



2 In this case, we need an additional relation v m j(tp)/v c j(tp) — ("fm/'fc) 2 oc (e e eB,r^t p ni) 2 for the 
magnetization estimate where 7 m and 7 C are the random Lorentz factors of electrons corresponding to the 
typical and cooling break frequencies, respectively. 
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4. Shocks in the intermediate regime 

The simple estimates of r d and Yd, which we have discussed in section |2j provide useful 
insights into the fireball dynamics. However, these are order-of-magnitude estimates, and 
obtained by assuming that the reverse shock is ultra-relativistic or Newtonian. Since most 
observed bursts are actually in the intermediate regime £o ~ 1> we here consider a better 
approximation which is similar to one discussed by Nakar & Piran (2004). 

The deceleration of an expanding shell happens when it gives a significant fraction of 
the kinetic energy to the ambient medium. Equalizing the energy in the shock ambient 
matter with E/2, we obtain = 2~ 1/ ' 3 //T 2 / 3 . The Lorentz factor Td in the shock regions is 
given as a function of the initial Lorentz factor Tq and the density ratio n^/ni between the 
unshocked shell material and the ambient medium (Sari & Piran 1995). For a homogeneous 
shell with width A, the particle density is ~ (E/m p c?T())/(47rR d AT()). Since the shock 
jump conditions and equality of pressure along the the contact discontinuity give a relation 
Hi/ Tii ~ ^T\l \(4r d + 3)(T d — 1)] , where T d , is the Lorentz factor of shocked shell material 
relative to the unshocked shell. We get an equation [f| for x = T d /T as 

24 T 8/3 

£ 2 (7) 

q 2 2 / 3 (l-x) 2 (2 + 3x + 2x 2 )' 1 ' 

where £ = (//A) 1 / 2 r o 4 ^ 3 and we have used ~ {x + l/x)/2. The corresponding results 
are shown in figure [2J For £ <C 1, we obtain x ~ 0.47£ 3//4 , while for £ ^> 1, we obtain 
x ~ 1. In the rest of the paper, we call the estimates obtained in the this section as the 
approximation (J7J) or the estimates based on Equation ([7]), while the estimates discussed 
in section 2 are called the conventional estimates. Since Equation (J7J) gives x for a given 
£, to estimate x at the deceleration radius r d , we need to use the value of £ at r d . In the 
thick shell regime £o < 1, £ is a constant during the deceleration process and we can use 
the initial value £o- However, in the thin shell regime £o 3> 1, due to the shell's spreading 
A ~ A + r d /Tl ~ A (1 + £q), the value £ ~ £ (l + £o) -1 ^ 2 is always about unity at r d (Sari & 
Piran 1995). Then, if we plot x and the reverse shock temperature (f d — 1) as functions of the 
initial value £o, they are expected to flatten in the intermediate regime. Since £ ~ £o in the 
thick shell and intermediate regime, we can directly compare the two approximations, and 
we find that the conventional approximation overestimates x and (f^ — 1) in the intermediate 
regime. The conventional estimates are x = $ 4 and (Td — 1) = £o~ 3 ^ 4 /2 for £ < 1 and x ~ 1 
and (f d - 1) ~ 1 for £ > 1. 

Using the deceleration radius r d and the Lorentz factor Td, the deceleration time is 



3 Assuming rd — Z 3/ ' 4 A 1 / 4 , Nakar & Piran (2004) have obtained a similar equation. 
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td ~ r rf/(2cr^) = l/(2^ 3 cT 8 d ). For the solution of Equation ([7]), we have an estimate of 
the correction factor Ct = 2~ 4 / 3 x _8//3 . Assuming no gradients in the distribution functions 
of the pressure and velocity in the shock regions, we obtain v m>r /v m ,f ~ (Xd — l) 2 / 1 ^ and 
F v ,max,r/ F Ujm ax, f ~ F 2 d /T where we have assumed R B = e B>r /e B j = 1- Then, we get the 
correction factors C m ~ (1 — x) 4 /(4x 4 ) and Cf ~ x 2 . 



5. Numerical Simulation 

The two analytic estimates which we have discussed include approximations (e.g. a 
simplified shock approximation and no gradients in the distribution functions of hydrody- 
namics quantities in the shock regions). Furthermore, the estimate (JTj) gives the Lorentz 
factor for a given £ at the deceleration radius r^, instead of the initial value £o- Since 
the typical frequency of the reverse shock emission is sensitive to the temperature (T^ — 1), 
it is important to investigate how £ at depends on £ ( or where £ becomes a constant) 
and what asymptotic value the reverse shock temperature takes in the thin shell regime. To 
examine the accuracy of the approximations and evaluate the shock Lorentz factors and the 
corrections factors C t , C m and Cf, we employ a spherical Lagrangian code based on the 
Godunov method with an exact Riemann solver (Kobayashi et al. 1999; Kobayashi & Sari 
2000; Kobayashi & Zhang 2007). No MHD effects are included in our purely hydrodynamic 
calculations. However, if the magnetization of a fireball is not too large (i.e. the ratio of 
magnetic to kinetic energy flux a ^ 0.1; Giannios et al. 2008; Mimica et al. 2009), the 
dynamics of shocks is not affected by magnetic fields, and our numerical results can be used 
to model the synchrotron emission from forward and reverse shocks. We will evaluate the 
correction factors for R B = e B ^/e B j = 1. 

The initial configuration for our simulations is a static uniform fireball surrounded by 
a uniform cold ISM. The hydrodynamic evolution is evaluated through the stages of initial 
acceleration, coasting, energy transfer to the ISM and deceleration. The evolution of a fireball 
is fully discussed in Kobayashi et al. (1999). We assume explosion energy E Q = 10 52 ergs and 
ambient density ri\ = 1 proton cm -3 for all the simulations, while we vary the dimensionless 
entropy (40 < Tq < 10 3 ) and the initial fireball size (ro = 10 9 cm, 6 x 10 11 cm, or 
6 x 10 12 cm) to cover a wide rage of £o- 
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5.1. Spectra and Light Curves 

We evaluate shock emission as a sum of photons from Lagrangian cells (fluid cells) 
in numerical calculations. First consider a single fluid cell with Lorentz factor T, internal 
energy density e, and particle density n and mass m in a shocked region (forward or reverse 
shock). Electrons are assumed to be accelerated in the shock to a power-law distribution 
with index p = 2.5 above a minimum Lorentz factor j m . We assume that constant fractions 
e e = 6 x 1CT 2 and e B = 6x 1CT 3 of the shock energy are given to electrons and magnetic fields, 
respectively. Our results are insensitive to the exact values of the microphysical parameters 
as long as Rb = 1, but they are included here for completeness. The typical random Lorentz 
factor and the energy of magnetic fields evolve as 7 m oc e/n and B 2 oc e. The typical 
synchrotron frequency in the observer frame is v m oc T^B, and the peak spectral power 
is F u , max oc N e TB for a total number N e ~ m/m, p of electrons in the cell. As we use a 
Lagrangian code, N e remains constant throughout the numerical evolution. The flux at 
a given frequency above v m is F v = F^ max (u /V™)"^ 1 ^ 2 oc mT^ p+1 ^ 2 e^ 5p ' 3 ^ 4 n 1 ~' p , while 
below v m we have a synchrotron low-energy tail as F u = F^ max (v /um) 1 ^ 3 oc mr 2 / 3 e _1 / 3 n 2 / 3 . 
Then, the emission from a cell can be estimated by using hydrodynamic quantities. 

We treat a fluid cell as a particle that continually emits photons. However, we only have 
the locations of the cell {rj} and flux estimates {F u j} at discrete timesteps {ti a b,j} where 
the subscript j indicates quantities at lab timestep j and rj indicates the inner boundary of 
the cell. We assume all the photons are emitted from the inner boundary (i.e. we neglect 
the radial width of the cell). Prior to the light curve construction we generate a series of 
(logarithmic) bins with boundaries {t^} in the observer time domain, and we assume bin k 
is bounded by tk and tk+i- We now consider the emission from a single fluid cell between 
two consecutive lab timesteps: j and j + 1. Since a photon emitted at timestep j arrives 
at the observer at tj = tiabj — rj/c, photons spread over observer time bins between tj and 
tj+i. Note that the observational time tj monotonically increases with j. Assuming that the 
observed flux F u evolves linearly between tj and and that the observer detects photons 
between v R and v R + dis, we can estimate the amount of energy deposited to each time bin. 

By monitoring the entropy evolution of a fluid cell, we can determine when the cell is 
heated by a forward or reverse shock. Then, we take into account all the timesteps after the 
shock heating for the construction of the light curve of the single fluid cell. We can apply this 
technique to all the cells inside (or outside) the contact discontinuity, and the total energy 
from all the cells in each time bin is divided by the bin size to get the reverse (or forward) 
shock light curve. It is then the simple matter of finding the maximum flux to obtain the 
peak time t p of the reverse shock emission. 

To numerically define a property / of the fireball shell at the peak time, we consider an 
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average value 

(/> = • (8) 

where the summation is taken over all the fluid cells {i} which are inside the contact discon- 
tinuity (i.e. in the reverse shock region) and which have contributed to the peak flux and 
SEi is the the contribution from fluid cell i to the peak time bin. 

At the peak time, we construct the spectrum. For this purpose, we set up a series of 
bins {v q } in the frequency domain. For the peak time bin (i.e. the time bin which gives 
the peak flux), we know which fluid cells have contributed, and at which lab timestep it 
has happened. Let us assume that a fluid cell deposits energy between lab timestep j and 
j + 1. Assuming a liner evolution of the luminosity J^ q+1 F v dv between the timesteps, we 
can estimate how much energy is deposited in each frequency bin at the peak time (the peak 
time bin). After summing up all the the energy deposited by the relevant fluid cells in each 
frequency bin, we divide the energy by the frequency bin size to get the spectrum at the 
peak time. 



5.2. Comparison of the Estimates and the Correction Factors 

Figure [2] shows the Lorentz factor Y & and the reverse shock temperature (f ^ — 1) at 
the peak time t p . For the numerical results (the dots), we have used Equation (JSJ) with 
(Td — 1) = e/nrripC 2 to obtain the average values, and we have assumed A = r to estimate 
£o- The numerical results and the conventional approximation are plotted against £o> while 
the approximation ([7]) is plotted against £. As we have discussed in Section 4, when the 
initial value is high ^> 1, the £ parameter is expected to decrease to order-of-unity during 
the evolution. One finds that such flattening in the numerical results occurs at a rather high 
value £ ~ several. The approximation ([7]), especially (1^—1) in the intermediate regime and 
in almost the whole range, is in better agreement with the numerical results, compared 
to the conventional estimates. The green dashed-dotted line in the bottom panel indicates 
the numerical asymptotic value: (Td — l) ~ 8 x 10~ 2 . 

Using a numerical peak time t p , we estimate the correction factor C t = ct p T^ 3 /I as 
a function of £o- The results are shown in the top panel of Figure [3J The conventional 
approximation well explains the numerical results in the thick shell regime £o < 1 but h 
breaks down in the thin shell regime £q > 1. In the thin shell regime, the numerical Ct is 
lower by a factor of ~ 5 than the conventional estimate which is equivalent to the numerical 
peak time being earlier than expected. Since for simplicity we have neglected a factor of 2 
in the conventional estimate as td = 1/cT^ 3 instead of l/2cT^ 3 , C t = 0.5 would be more 
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appropriate for the conventional estimate in the thin shell regime. However, the numerical 
results are still smaller. This is in part due to the gradients in the distribution functions 
of hydrodynamics quantities in the reverse shock region. The numerical distributions have 
a higher value at the contact discontinuity, and it decreases toward the tail (see Figure 3 
in Kobayashi & Zhang 2007). It makes the contribution of photons from the inner parts of 
the fireball shell less significant, reducing the effective width of the emission region in the 
shell, and the shock emission peaks earlier than in the homogeneous shell case. Since as we 
will see later, the magnetization estimate is rather insensitive to C t (and the peak time), we 
discuss only the line-of-sight emission in this paper. However, it is possible to include the 
high latitude emission at expense of computational time, and we have obtained very similar 
results for several selected cases. With the addition of the high latitude emission the overall 
light curve appears smoother with slightly shallower decay features. The position of the peak 
time increases by ~ 50%. 

Using the numerical values of the typical frequency ratio at the peak time, we estimate 
the correction factor C m = Tq [u m . r {tp) / u m j{t p )} as a function of £o- The results are shown 
in the middle panel of Figure [3J The conventional estimate is in good agreement with 
the numerical results in the thick shell regime, but as we expect from C m oc {T d — l) 2 , it 
overestimates C m by a factor of ~ 10 2 in the thin shell regime. Finally the bottom panel 
of Figure [3] shows the results for Cp = r^ 1 [F^ max>r {tp) / F u>max j{t p )]. The conventional 
approximation overestimates the amount of flux emitted by the reverse shock especially in 
the intermediate regime as Nakar & Piran (2004) have pointed out. The estimates based 
on Equation (J7J) provide a better approximation for all the three correction factors in the 
intermediate regime. The red dashed lines in the three panels indicate the numerical fitting 
formulae C t = N t + £ ~ 2 with N t ~ 0.2, C m = N m + ^ 3 with iV m ~5x 10~ 3 and Cp 1 = 
Np + M f £ Pf with Np ~ 1.5, Mp ~ 5 and P F ~ 1.3. 

Figure S] illustrates wide band spectra at the peak time. We here consider three numer- 
ical cases with £o — 0.1, 1 or 10. The black line indicates the conventional estimate in which 
the typical frequency v m ^ r of the reverse shock emission is lower by a factor of Y\ than that of 
the forward shock emission, and the peak flux F v ^ maxr of the reverse shock emission is higher 
by a factor of T . However, our numerical results show that in the thin shell regime v m>r is 
lower further by a factor of ~ 10 2 than the conventional estimate (the red line), and that 
in the intermediate regime the peak flux F Vimax , r would be lower by a factor of several (the 
green line). These indicate that the reverse shock emission would be elusive if the typical 
frequency of the forward shock is at around the optical band and if the forward and reverse 
shock have the same microscopic parameters (Nakar & Piran 2004; Mimica et al. 2010; 
Melandri et al. 2010). In the thick shell regime, the peak frequency of the reverse shock 
emission is closer to that of the forward shock emission u m>r /v m j ~ £^ 3 rQ 2 (the blue line). 
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We might have a better chance to detect the reverse shock component in early afterglow, 
although the light curve peaks earlier than in the thin shell regime. 



5.3. Initial Lorentz Factor and Magnetization Parameter 

The initial Lorentz factor can be evaluated by using the peak time t p ~ t^, 



C t l 



3/8 



where I is basically a known quantity from the prompt gamma-ray and late-time afterglow 
observations, the estimate depends very weakly on more fundamental parameters T oc 
/ 3//§ oc (£ , /n 1 ) 1 / 8 , and we had obtained numerical result C t ~ 0.2 + In principle, we 
can estimate £o from observable quantities. Since the duration T of the prompt gamma-rays 
gives a rough estimate of the width A ~ cT (Kobayashi et al. 1997), using Equation ([9]) and 
the numerical C f (£ ), we obtain £ ~ 5 1 / 2 a/ \t p /T) — 1 and C t ~ 0.2(1 — T /t p )~ l . In the thin 
shell regime, an early afterglow peaks well after the prompt gamma-ray emission (Sari 1997), 
and we have C t ~ 0.2. However, in the thick shell regime, the peak time is almost equal to 
the width Aq/c. The approximation Ao = cT might not be accurate enough to discuss the 
exact value of £o- Since the flux before the optical peak t p is sensitive to the initial profile of 
the shell and in particular to the rising index of the light curve might be used to break 
the degeneracy of the £o estimate in the thick shell regime. Nakar & Piran (2004) have 
numerically estimated the rising index for a homogeneous shell in a rage of 0.05 < £o < 5 
as a rise ~ 0.6 [1 + p(£ — 0.07£ 2 )]. A slow (rapid) rise is a signature of the thick (thin) shell 
regime (Kobayashi 2000). 

The magnetization parameter Rb can be estimated by using T R t , and Rp. For the 
typical decay index of the reverse shock emission a = 2, the conventional estimate is Rb,cou = 
(R 3 F T /R3) 2 / 5 in the thin shell regime where To = (l/ct p ) 3 ^ 8 . Then, we obtain a correction 
factor for the magnetization parameter as 

D ^3/20 

Rb Ct (10) 



E> ^6/5^4/5 



One finds that the estimate is rather insensitive to C t . In Figure [51 the numerical results 
are plotted with the approximations. For a typical GRB (£q ~ 1), the conventional approx- 
imation (black solid line) underestimates the magnetization parameter by a factor of ~ 10. 
A more extreme discrepancy occurs in the thin shell regime, and the magnetization param- 
eter is underestimated by a factor of ~ 10 2 . The estimate based on Equation (J7|) (black 
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dashed line) describes the numerical results reasonably well in the intermediate and thick 
shell regime. 

6. Case Studies 

To demonstrate our scheme, we discuss early optical afterglows detected from GRB 
990123 and GRB 090102. 

GRB 990123: This burst is one of the brightest GRBs observed so far. The basic 
parameters include z = 1.6, E ~ 1.4 x 10 54 ergs, and T ~ 60s. (e.g. Kobayashi & Sari 2000 
and references therein). The gamma-ray profile is dominated by two pulses, each lasting ~ 8s, 
separated by 12s. A bright optical flash is detected during the prompt emission (Akerlof et 
al. 1999), the optical emission peaks at t p ~ 50s with 9 mag, and rapidly decaying with a 
power-law slope of a ~ 2 and it becomes shallower at a late time t p j < 0.1 days. We have 
only one optical data point before the peak, and it provides a lower limit of the rising index 
a r i se > 2 and the the peak flux. We conservatively assume that the peak flux is 9 mag. Since 
the optical peak is comparable to the GRB duration (especially to the duration of the main 
two pulses) and the rising is rapid, this is an intermediate case £o ~ 1, with the correction 
factors C t ~ 1.2, C m ~ 1 and C F ~ 0.16. The initial Lorentz factor is about T ~ 460%^ 8 . 
Assuming t p j ~ O.ldays, one has R t ~ 170 and R F ~ 5000, therefore the magnetization 
parameter is R B ~ 6000. Since Zhang et al. (2003) found a magnetization parameter of 
Rb ~ 15 2 = 225 based on the conventional approximation with £o ~ 1 5 the magnetization 
estimate with the corrections is larger by an order of magnitude. 

GRB 090102: This burst shows a significant polarization at the 10% level in the 
early optical afterglow (Steele et al. 2009). The basic parameters include z = 1.5, and 
E ~ 5.8 x 10 53 ergs (e.g. Gendre et al. 2010 and reference therein). The prompt gamma-ray 
emission lasts for 27s and comprises four overlapping peaks starting 14s before the GRB 
trigger. The optical light curve, beginning with 13 mag at 40s after the GRB trigger, can 
be fitted by a broken power-law whose flux decays as a function of time (F oc t~ a ) with a 
gradient a ~ 1.5 that then flattens to a ~ 1 after ~ 10 3 s. If we assume that the optical 
emission peaks at the first data point (the mid time is ~ 60s after the beginning of the GRB) 
and t p j ~ 10 3 s, we obtain £ ~ 2.4 and the correction factors: C t ~ 0.3, C m ~ 8 x 10~ 2 , and 
C F ~ 0.3. Using R t ~ 17, R F ~ 70 and T ~ 240n^ 1/s , we obtain R B ~ 97. On the other 
hand, if £o — 0.1, we obtain T = 2200n ] ~ 1 ^ 8 and R B ~ 42. The detection of the significant 
polarization in the steep decay phase (t = 160 ~ 220s) and the large magnetization parameter 
Rb indicate the presence of large-scale magnetic fields in the fireball. 
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7. Conclusions 

We have discussed a revised method to estimate the magnetization degree of a GRB 
fireball. We use the ratios of observed properties of early afterglow so the poorly known 
parameters for the shock microphysics (e.g. e e and p) would cancel out. Since the estimate 
depends only weakly on the explosion energy and the fireball declaration time, the estimate 
does not require the exact distance (redshift) to the source as an input parameter. Since most 
observed events fall in the intermediate regime between the thin and thick shell extremes, 
we have provided a new approximation for the spectral properties of the forward shock and 
reverse shock emission, which well describes the numerical results in the intermediate regime. 
The previous standard approach underestimates the degree of fireball magnetization by a 
factor of 10 ~ 100. 

In the GRB phenomena, extreme relativistic motion with To > 100 is necessary to 
avoid the attenuation of hard gamma-rays. The acceleration process is likely to induce a 
small velocity dispersion inside the outflow Ar ~ T (e.g. thermal acceleration). If internal 
shocks are responsible for the production of the prompt gamma-rays, the dispersion should 
be even larger (Beloborodov 2000; Kobayashi & Sari 2001). The velocity dispersion leads 
to the spreading of the shell structure in the coasting phase and the £ parameter decreases 
as £ oc A -1 / 2 . As an order-of-magnitude estimate, when the initial value £ > 1, the reverse 
shock always becomes mildly relativistic (£ ~ 1) at the declaration radius and the reverse 
shock temperature (f d — 1) is expected to be insensitive to the initial value £o- However, 
it is difficult to analytically quantify the asymptotic reverse shock temperature. We have 
numerically shown that the spreading effect becomes significant at rather high values £o £ 
several, and that the asymptotic value is (T d — 1) ~ 8 x 10~ 2 . 

We have confirmed that, especially in the intermediate regime £o ~ 1, the reverse shock 
emission is much weaker than the standard estimates, and that in the thin shell regime the 
typical frequency of the reverse shock emission is much lower than the standard estimates 
(Nakar & Piran 2004). If the fireball is not magnetized Rb = £B,r/ e B,f — 1, the reverse 
shock emission more easily falls below the forward shock emission. The lack of optical 
flashes from most GRBs might be partially explained in a revised non-magnetized model. If 
the fireball shell does not spread even in the thin shell regime (i.e. the velocity distribution 
is completely uniform), only a small fraction of the kinetic energy of the shell is converted to 
thermal energy in the reverse shock, and the reverse shock emission is practically suppressed 
in the thin shell regime £o ^ 1- 

The broadband afterglow emission of GRB 990123 is modeled to find €bj ~ 4 x 10~ 7 
(Panaitescu & Kumar 2004). Although there are no estimates available for GRB 090102, the 
broadband modeling generally shows that it is in a range of e B j = 10~ 5 — 10 _1 (Panaitescu 
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& Kumar 2002). Using the estimated values for the magnetization parameter this would 
put the reverse shock magnetic energy fraction as es, r ~ 2 x 10~ 3 for GRB 990123 and 
€B, r ~ — 10 for GRB 090102 with £ = 2.4. Although magnetic pressure would suppress 
the formation of a reverse shock if ee^^O.l (Giannios et al. 2008; Mimica et al. 2009), 
the low values are consistent (or the result includes the parameter region consistent) with 
the basic assumption in our analysis (i.e. magnetic fields do not affect the reverse shock 
dynamics). If a future event indicates a high value e£, r >0.1, an interesting possibility to 
reconcile the problem is that the prompt optical emission (and prompt gamma-rays) would 
be produced through a dissipative MHD processes rather than shocks (Giannios & Spruit 
2006; Lyutikov 2006; Giannios 2008). 

Our magnetization estimates are slightly lowered if the blast wave radiates away a 
significant fraction of the energy in the early afterglow E oc t~ 17 ^/ 16 (Sari 1997). For 
e e = 0.1, the blast wave energy becomes smaller by a factor of 1.7 between t = 50s and 
O.ldays, by a factor of 1.3 between t = 60s and 10 3 s, then the estimates of the magnetization 
parameter are reduced by a factor ~ 2 for GRB 990123 and ~ 1.5 for GRB 090102. These 
afterglows show a flattening after the steep decay phase. Since the reverse shock emission 
dominates at early times, the forward shock peak time is rather uncertain in the both cases. 
Our magnetization estimates are lowered if the forward shock peaks much earlier than the 
flattening time. As we have discussed in section 3, the estimate depends on the peak time 
as R B oc t° p f for GRB 990123 and t° p f for GRB 090102. If we assume that the forward 
shock emission peaks at the deceleration time t p j ~ 50s for for GRB 990123 and ~ 60s for 
GRB 090102, the magnetization parameter would be R B ~ 3200 and ~ 12, respectively. 
The results still suggest magnetized fireballs. 

We thank Ehud Nakar and Elena M. Rossi for useful discussion. This research was 
supported by STFC fellowship and grant. 
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Fig. I. — Early steep decay: optical afterglow, produced as a composition of the reverse 
shock emission (solid line) and forward shock emission (dashed line). Two peaks (top panel) 
and a flattening (bottom panel) in the light curve. 
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Fig. 2. — r^/To (Top panel) and (f^ — 1) (Bottom panel) as functions of £o or £. The 
conventional estimates (black solid lines), the approximation (J7J) (black dashed lines), and 
the numerical results (orange dots: r = 10 9 cm, blue dots: 6 x 10 11 cm, and red dots: 
6 x 10 12 cm). See the text concerning to the choice of £o or £. 
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Fig. 3. — The correction factors as a function of £o or £■ Top panel: The normalized peak 
time C t , Middle panel: The normalized typical frequency ratio C m , and Bottom panel: The 
normalized peak flux ratio Cf- The conventional estimates (black solid lines), the estimate 
based on the approximation (J7J) (black dashed lines), the numerical results (blue, red and 
oranee dots are the same as in Fieurel2l), and numerical fitting formulae (red dashed lines"). 




Fig. 4. — Numerical wide band spectra at the peak time: £o = 0.1 (blue line), 1 (green 
line) and 10 (red line). The black line shows the conventional estimate (thin shell case). 
The frequency and flux are normalized by the typical frequency and peak flux of the forward 
shock emission, respectively. T = 300 and A ~ 3 x 10 9 cm (red line), ~ 3 x 10 11 cm (green 
line), or ~ 3 x 10 13 cm (blue line) are assumed. 
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Fig. 5. — The correction factor for the magnetization parameter RB/RB,con as a function 
of £ or £• The legend is the same as figure [3] with the best fit equation (red dashed line) 
obtained by the combination of the best fits to C t , C m and Cp. a = 2 is assumed. 



